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We predict that unpolarized charge current driven through the longitudinal leads attached to 
ballistic quantum-coherent two-dimensional electron gas (2DEG) in semiconductor heterostructure 
will induce a pure spin current, which is not accompanied by any net charge flow, in the transverse 
voltage probes. Its magnitude can be tuned by the Rashba spin-orbit (SO) interaction and, more- 
over, it is resilient to weak spin-independent scattering off impurities within the metallic diffusive 
regime. While the polarization vector of the spin transported through the transverse leads is not 
orthogonal to the plane of 2DEG, we demonstrate that only two components (out-of-plane and lon- 
gitudinal) of the transverse spin current are signatures of the spin Hall effect in four-probe Rashba 
spin-split semiconductor nanostructures. The linear response spin Hall current, obtained from the 
multiprobe Landauer-Biittiker scattering formalism generalized for quantum transport of spin, is 
the Fermi-surface determined nonequilibrium quantity whose scaling with the 2DEG size L reveals 
the importance of processes occurring on the spin precession mesoscale Lso (on which spin precesses 
by an angle tt) — the out-of-plane component of the transverse spin current exhibits quasioscillatory 
behavior for L < Lso (attaining the maximum value in 2DEGs of the size Lso x Lso), while it 
reaches the asymptotic value in the macroscopic regime L 3> Lso- Furthermore, these values of the 
spin Hall current can be manipulated by the measuring geometry defined by the attached leads. 

PACS numbers: 72.25.Dc, 73.23.-b, 85.75.Nn 



I. INTRODUCTION 

Current efforts in spintronics are to a large extent di- 
rected toward gaining control of electron spin in semicon- 
ductor structures, which are ubiquitous in conventional 
electronics, and exploiting it as a carrier of classical or 
quantum information.^ Although spin physics in semi- 
conductors is an old subject (2*^ spintronics has reignited 
interest in the role of spin-orbit (SO) couplings in con- 
densed matter systems. While they originate from rela- 
tivistic corrections to the Schrodinger equation, SO in- 
teractions for itinerant electrons in semiconductor nanos- 
tructures can be much stronger than for particles mov- 
ing through electric fields in vacuumA^ Therefore, they 
are envisaged to be a tool for all-electrical^ spin current 
generation and manipulation, where electric fields can be 
produced to control electron spin in far smaller volumes 
and on far shorter time scales than it is possible with 
conventional magnetic field-based spin controli 

Recent theoretically unearthed intrinsic spin Hall ef- 
fect in hole doped^ (such as bulk p-GaAs or p-Ge de- 
scribed by the Luttinger effective Hamiltonian for heavy 
and light holes), or electron dopedA (such as 2DEG in 
n-type heterostructures with structure inversion asym- 
metry which gives rise to the Rashba-type of SO cou- 
pling^'"''^) semiconductors suggests that pure (i.e., not ac- 
companied by any net dissipative charge current) trans- 
verse spin current could be generated in these systems 
due to longitudinal electric field. The correlation be- 
tween spin orientation and carrier velocity, induced in 
this effect in the presence of an external electric field, 
essentially requires some type of SO coupling which is 
strong enough to spin-split the Bloch energy bandsi 

On the other hand, it has been known for a long time^'^ 




FIG. 1: (Color online) The four-probe mesoscopic bridge for 
the detection of pure spin Hall currents. The central region 
is 2DEG where electrons are confined within a semiconductor 
heterostructure by an electric field along the 2-axis which in- 
duces the Rashba SO coupling. The four attached leads are 
clean, non-magnetic, and without any SO coupling. The un- 
polarized (7i = 0) charge current (Ji 7^ 0) through the longi- 
tudinal leads induces spin Hall current in the transverse leads 
which act as the voltage probes V2 = V3 0, I2 = I3 = 0. We 
predict that the polarization vector of the spin transported by 
pure {I2 = -^^2 + -^^2 — 0) spin current /| = ^(-^2 ~ ^2) is not 
orthogonal to the plane of 2DEG. 



that SO dependent scattering off impurities will deflect 
spin-t (spin-|) electrons predominantly to the right (left), 
thereby giving rise to the extrinsic spin Hall effect where 
pure spin current flows perpendicular to the longitudi- 
nal unpolarized charge transport. However, the intrin- 
sic spin Hall current is expected to be several order of 
magnitude larger than the extrinsic one, thereby promis- 
ing all-electrical solution to spin injection problems* and 
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opening new avenues for semiconductor spintronics ap- 
plications,''^ 

While the properties of the intrinsic spin Hall current 
have been delineated through the semiclassical analysis of 
infinite homogeneous SO coupled semiconductor systems 
in the clean limit guiding experimental detection of 
such effects requires a quantitative prediction for the spin 
current flowing through the leads attached to a finite-size 
sample, as exemplified by the bridge in Fig. ^ This is 
analogous to profound developments in our understand- 
ing of quantum Hall effect ensuing from the comparison 
of the macroscopic charge transport in bulk samples with 
the gap in the energy spectrum to the mesoscopic trans- 
port through gapless chiral edge states of multiterminal 
bridges employed in experiments; -'-^ 

For example, within a finite-width strip no charge 
or spin current can flow across its boundaries, so that 
nonequilibrium spin accumulation^^ will appear near the 
lateral edges to generate compensating current in the di- 
rection opposite to the spin Hall current. Very recent ex- 
perimentsi^ have indeed demonstrated for the first time 
the existence of such spin accumulation, which has oppo- 
site sign on the lateral edges of the wire, as the manifes- 
tation of the spin Hall effect(s) in two-terminal devices. 
Thus, when ideal (i.e., spin and charge interaction free) 
transverse leads are attached at the lateral edges of the 
2DEG region in Fig. ^ pure {I^ + I2 =0) spin current 
^2 = 55(^2 ^ ^2) should emerge in the probe 2 of the 
bridge. 

Here we predict a novel type of effect which exhibits the 
spin Hall phenomenology in mesoscopic finite-size struc- 
tures. The transverse pure spin current /| in Fig. ^ is 
induced by injecting unpolarized charge current through 
the longitudinal leads into the SO coupled central region 
which contains no impurities. The central region is as- 
sumed here to be a finite-size Rashba spin-split 2DEG, 
which is phase- coherent (i.e., electron is described by a 
single wave function within the sample) and in the ballis- 
tic transport regime where electrons do not feel any elec- 
tric field while propagating through the 2DEG. Within 
the 2DEG carriers are subjected to the Rashba type of 
SO coupling, which is described by the following effective 
mass Hamiltonian 

^ 2 

i? = ^ + |('3-xp).z + l4onf(y), (1) 

that takes into account structure inversion asymmetry2i4 
(of the confining electric potential and differing band dis- 
continuities at the heterostructure quantum well inter- 
face^). Here p is the momentum operator in 2D space, 

= (o'a;, CTy, o'z) is the vcctor of the Pauli spin matrices, 
a is the strength of the Rashba SO coupling;2i^ and 
V^onf (2/) is the transverse confining potential. 

The paper is organized as follows. In Sec. we recast 
the Landauer-Biittiker multiprobe charge current formu- 
las in terms of spin-resolved charge current, which then 
allows us to introduce the multiprobe spin current for- 
malism that yields the general expression for the linear 



response spin Hall conductance of the bridge in Fig. ^ 



in terms of the spin-resolved transmission probabilities 
between different leads. We apply this formalism to 
a perfectly clean bridge to demonstrate in Sec. IIIII the 
existence of three non-zero spin conductances (corre- 
sponding to three components [/|]^, [li]^ , [liY of the 
vector of pure spin current flowing through the trans- 
verse leads), two of which GJ^ = [/2]^/(^i ~ ^) ^^id 
G^fj — [l2]^/{Vi — V4) represent the signature of the 
mesoscopic spin Hall effect in Rashba spin-split struc- 
tures. They are determined by the density of electrons 
and the Rashba SO coupling strength a, as well as by the 
measuring geometry of the whole device (i.e., interfaces, 
boundaries, and the attached electrodes). Furthermore, 
we find two different scaling laws for and GJ^, de- 
pending on whether the device is smaller or greater than 
the mesoscale set by the spin precession length Lso on 
which spin precesses by an angle tt. Section IIVI shows 
that the predicted effect is able to survive weak disor- 
der (introduced as spin-independent scattering off static 
impurities) — the spin Hall conductances G^^ and Gf^ 
gradually diminish from their maximum values (set in the 
clean limit of Sec. lHIl) within the metallic diffusive regime 
and become negligible as the disorder is increased, but 
before the onset of strong localization effects is reached 
in phase-coherent SO coupled 2D structures. 

The magnitude of the intrinsic spin Hall effect in infi- 
nite homogeneous systems is captured by the spin Hall 
conductivity a sH = jyl^x which relates pure spin current 
density jy, flowing in the y-direction and carrying spins 
polarized solely along the z-axis, as a response to the lon- 
gitudinal externally applied electric field Ex- Using our 
results from Sec. IIHI and Sec. IIVI we analyze in Sec. IVlif 
any quantitative connection can be established between 
the bulk spin Hall conductivity <Jsh and the spin Hall 
conductance G^^u (i.e., between the intrinsic spin Hall 
current density jy, which is not conserved in the bulk 
and depends on the whole SO coupled Fermi sea, and 
[I2Y which we find to be a Fermi-surface quantity and 
conserved total spin current throughout the ideal leads). 
This analysis reveals different origins of the mesoscopic 
spin Hall effect, which is governed by the processes on 
the mesoscalai'^ Lso f^nd represents the nonequilibrium 
manifestation of SO couplings in confined ballistic semi- 
conductor nanostructures. We conclude in Sec. IVII 

II. SCATTERING APPROACH TO QUANTUM 
TRANSPORT OF SPIN CURRENTS IN 
MULTIPROBE GEOMETRIES 

The mesoscopic experiments on quantum Hall bridges 
in the early 1980s were posing a challenge for theoret- 
ical interpretation of multiterminal transport measure- 
ments.^" By viewing the current and voltage probes on 
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equal footing, Biittiker— has provided an elegant solu- 
tion to this problem in the form of a multiprobc for- 
mulcOiii^ 

which relates, via the conductance coefficients Gpq, 
charge current Ip = /J + in probe p to the volt- 
ages Vq in all other probes attached to the sample. To 
study the spin-resolved charge currents (cr =tj-l) of 
individual spin species "f, | we imagine that each non- 
magnetic lead in Fig. ^ consists of the two leads allow- 
ing only one spin species to propagate (as realized by, 
e.g., half-metallic ferromagnetic leads). Upon replace- 
ment Ip — !■ Ip and Gpq — > Gpq , this viewpoint allows us 
to extract the multiprobe formulas for the spin-resolved 
charge currents^^'"'^^ , thereby obtaining the linear re- 
sponse relation for spin current Ip — ^{Ip — Ip) flowing 
through the lead p 

~ 2^ ^ y K^qj) Gqp ~ Gqp ~ '-^gp)^ 
1 

-{Gll + Gll^-Gll-G-l^q)Vq]. (4) 

Below we simplify the notation by introducing the labels 
Qin ^ /^TT _|_ qU _ nil _ nil and G°"* = G^T _|_ nil _ 

pq PQ PQ pq ^ pq ^^^^ ^ pq ^ pq ' ^ pq 

Gj,^ — G^}j. Furthermore, these coefficients have transpar- 
ent physical interpretation: ■^G°p*'Vp is the spin current 
flowing from the lead p with voltage Vp into other leads q 
whose voltages are Vq, while -^Gp^Vq is the spin current 
flowing from the leads q ^ p into the lead p. 

The standard charge conductance coeBicientai^'i^^ in 
the multiprobe Landauer-Biittiker formalism Eq. ((3J) are 
expressed in terms of the spin-resolved conductances as 
Gpq = Gil + Gil + Gpq + Gl\. Their introduction in 
1980s was prompted by the need to describe linear trans- 
port properties of a single sample, with specific impurity 
arrangements and attached to specific probe configura- 
tion, by using measurable quantities (instead of the bulk 
conductivity which is inapplicable to mesoscopic conduc- 
torsi^). They describe total charge current flowing in 
and out of the system in response to voltages applied at 
its boundaries. 

Regardless of the detailed microscopic physics of trans- 
port, conductance coefficients must satisfy the sum rule 
Gqp — J2q Gpq iu Order to ensure the second equal- 
ity in Eq. ||2Jl, i.e., the charge current must be zero 
Vq = const. =^ Ip = m equilibrium. On the other hand, 
the multiprobe spin current formulas Eq. Q apparently 
posses a nontrivial equilibrium solution Vq = const. 
Ip ^ (found in Ref. ITt!) that would be entirely alien 
to the Landauer-Biittiker paradigm demanding usage of 
only measurable quantities. However, when all leads are 
at the same potential, a purely equilibrium non-zero term 
JsiG°pfVp - Gl^pVp) = !^{Gll - Gll)Vp becomes relevant 
for Ip [note that for devices in nonequilibrium the sum- 
mation in Eq. Q goes only over q ^ p leads] , canceling 



all other terms in Eq. to ensure that no unphysical 
total spin current Ip ^ can appear in the leads of an 
unbiased (V;j=const.) multiterminal devicei ^^i^" 

At zero temperature, the spin-resolved conductance co- 
efficients Gpq = X Sij |ti7crcr' where summation is 
over the conducting channels in the leads, are obtained 
from the Landauer-type formula as the probability for 
spin-cr' electron incident in lead q to be transmitted to 
lead p as spin-cr electron. The quantum- mechanical prob- 
ability amplitude for this processes is given by the matrix 
elements of the transmission matrix t''', which is deter- 
mined only by the wave functions (or Green functions) 
at the Fermi energy.^^ The stationary states of the struc- 
ture 2DEG + two leads supporting one or two conduct- 
ing channels can be found exactly by matching the wave 
functions in the leads to the eigenstates of the Hamilto- 
nian Eq. thereby allowing one to obtain the charge 
conductance from the Landauer transmission formula.— 
However, modeling of the full bridge geometry with two 
extra leads attached in the transverse direction, as well 
as existence of many open conducting channels, requires 
to switch from wave functions to some type of Green 
function formalism. 

For this purpose we represent the Rashba Hamilto- 
niani^ of the 2DEG in Fig. ^ in a local orbital basis 
defined on the L x L lattice (with lattice spacing a) asS^ 

H= K]e„,|m)(m| -io E |m)(m'| 

\ ^ (m.m') J 

+ ^{Py<^^x-Px<^^y)- (5) 

Here to (a unit of energy) is the nearest-neighbor hop- 
ping between s-orbitals (r|m) — ip{Y — m) on adjacent 
atoms located at sites m — {nix, my) of the lattice. 
Since momentum operator in the tight-binding represen- 
tation is (m|pa;|m') = 5jn'^^rn^±iih{m3: - m'^) /2a'^, the 
Rashba SO term (in which (g) stands for the tensor prod- 
uct of operators) introduces the SO hopping energy scale 
^so = ct/2a. In a perfectly clean 2DEG the on-site poten- 
tial energy is Em — 0, while disordered 2DEG can be sim- 
ulated via a random variable Em £ [—W/2,W/2] model- 
ing short-range isotropic scattering off spin-independent 
impurities. 

For non-interacting particle which propagates through 
a finite-size sample of arbitrary shape, the transmission 
matrices 

= ^/-Im ±p<»is- Glq ■ -yZ-Im ® 

imSp === (6) 

between different leads can be evaluated in a numeri- 
cally exact fashion using the real(g)spin-space Green func- 
tions.^^ This requires to compute a single object, the re- 
tarded Green operator 
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which becomes a matrix (i.e., the Green function) when 
represented in a basis |m) (g) |cr) S H.o^'Hs introduced by 
the Hamiltonian Eq. Here \a) are the eigenstates of 
the spin operator for the chosen spin quantization axis. 
The matrix elements G''(m'cr'; m, cr) — (m', o''|G""|m, ct) 
yield the probability amplitude for an electron to prop- 
agate between two arbitrary locations m and m' (with 
or without flipping its spin a during the motion) inside 
an open conductor in the absence of inelastic processes. 
Its submatrix G^^, which is required in Eq. consists 
of those matrix elements which connects the layer of the 
sample attached to the lead p to the layer of the sam- 
ple attached to the lead q. The unit operators lo and 
Is act in the orbital Tio and the spin Hilbert spaces Tis, 
respectively, which comprise the Hilbert space of a single 
spinfuU particle Tio ® Ti-s (via tensor product of vector 
spaces). The self-energy ® (r-retarded, a- 

advanced, T,p = [Sp]^) account for the "interaction" of 
an open system with the attached four ideal semi-infinite 
leadsi^ p. 

A direct correspondence between the continuous effec- 
tive Rashba Hamiltonian Eq. [with parabolic energy- 
momentum dispersion] and its lattice version Eq. Q 
[with tight-binding dispersion] is established by selecting 
Ep close to the bottom of the band (where tight-binding 
dispersion reduces to the quadratic one), and by using 
to = S^/(2m*a^) for the orbital hopping which yields the 
effective mass m* in the continuum limit. We elucidate 
further the connection between the standard effective 
Rashba Hamiltonian in continuous representation Eq. 
and its lattice version by interpreting the tight-binding 
parameters in Eq. (jSJ for particular experimental realiza- 
tion of a 2DEG in semiconductor heterostructures. For 
example, the InGaAs/InAlAs heterostructure employed 
in experiments of Ref. is characterized by the effective 
mass TO* = 0.05too (toq is the free electron mass) and 
the width of the conduction band A = 0.9 eV, which sets 
to — A/8 = 0.112 meV for the orbital hopping parameter 
on a square lattice (with four nearest neighbors of each 
site) and a ~ 2.6 nm for its lattice spacing. Thus, the 
Rashba SO coupling of 2DEG formed in this heterostruc- 
ture, tuned to a maximum value^ a = 0.93-10"^"'^ eVm by 
the gate voltage covering the 2DEG, corresponds to the 
SO hopping tso/to — 0.016 in the lattice Hamiltonian 
Eq. ©. 



A. General expression for the spin Hall 
conductance 

Since the total charge current Ip = I^+lj, depends only 
on the voltage difference between the leads in Fig. ^ we 
set one of them to zero (e.g., V4 = is chosen as the ref- 
erence potential) and apply voltage Vi to the structure. 
Imposing the requirement I2 = h — for the voltage 
probes 2 and 3 allows us to get the voltages V2/V1 and 
V3/V1 by inverting the multiprobe charge current formu- 
las Eq. Finally, by solving Eq. Q for /| we obtain 



the most general expression for the spin Hall conductance 
defined by Eq. (0) 
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(8) 

This quantity is measured in the units of the spin conduc- 
tance quantum e/47r (as the largest possible GsH in the 
transverse leads with only one open conducting channel), 
which is the counterpart of a familiar charge conductance 
quantum e^//i (as the natural unit for spin-resolved con- 
ductance coefficients G^^ ). 

In contrast to the charge current which is a scalar quan- 
tity, spin current has three components because of the 
vector nature of spin (i.e., different "directions" of spin 
correspond to different quantum mechanical superposi- 
tions of It) and II) states). Therefore, we can expect that, 
in general, the detection of spin transported through the 
transverse leads of mesoscopic devices will find its ex- 
pectation values to be non-zero for all three axes. We 
indeed find in Sec. Illll and Sec. llVl that all three compo- 
nents of the spin current in the transverse leads 2 and 3 
are non-zero. However, their flow properties 
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(9b) 
(9c) 



show that only the z- and the x-components represent 
the spin Hall response for the Rashba SO coupled four- 
terminal bridges. That is, if we connect the transverse 
leads 2 and 3 to each other (thereby connecting the lat- 
eral edges of 2DEG by a wire), only the spin current 
carrying z- and x-polarized spins will flow through them, 
as expected from the general spin Hall phenomenology 
where nonequilibrium spin Hall accumulation detected 
in experiments^^ has opposite sign^ on the lateral edges 
of 2DEG. 

Therefore, to quantify all non-zero components of the 
vector of transverse spin current in the linear response 
regime, we introduce three spin conductances G^j^ = 
mr/Vi, Glp = [/|]V^i, and GIh = [/|]V^i (assum- 
ing Vi — 0). They can be evaluated using the same 
general formula Eq. ((SJ where the spin quantization axis 
for t, I in spin-resolved charge conductance coefficients 
is chosen to be the x-, y-, or z-axis, respectively. For 
example, selecting (Jz|T) = +IT) and (Tz\i) = — |i) for the 
basis in which the Green operator Eq. |(7J is represented 
allows one to compute the z-component of the spin cur- 
rent [IpY- In accord with their origin revealed by Eq. lO, 
we denote G^jj and Gfj^ as the spin Hall conductances, 
while G^p is labeled as the "spin polarization" conduc- 
tance since it stems from the polarization of 2DEG by 
the flow of unpolarized charge current in the presence of 
SO coupling s^ ■ (see also Sec. IIII|I . 
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B. Symmetry properties of spin conductances 

Symmetry properties of the conductance coefficients 
with respect to the reversal of a bias voltage or the direc- 
tion of an external magnetic field play an essential role 
in our understanding of linear response electron trans- 
port in macroscopic and mesoscopic conductors^^*i^ For 
example, in the absence of magnetic field they satisfy 
Gpq = Gqp (which can be proved assuming a particular 
model for charge transport^^). Moreover, since the effec- 
tive magnetic field B/j(p) of the Rashba SO coupling de- 
pends on momentum, it does not break the time-reversal 
invariance which imposes the following property on the 
spin-resolved conductance coefficients G^^ = G~^ ~'' in 
the multiterminal SO coupled bridges>i2ii& 

In addition, the ballistic four-terminal bridge in Fig.Q] 
with no impurities posseses various geometrical symme- 
tries. It is invariant under rotations and refiections that 
interchange the leads, such as: (i) rotation C4 (C2) by an 
angle 7r/2 (tt) around the z-axis for a square (rectangular) 
2DEG central region; (ii) refiection a^x in the xz-plane; 
and (iii) reflection a^y in the j/z-plane. These geometrical 
symmetries, together with Gpq = Gqp property, specify 
V2/V1 = V3/V1 = 0.5 solution for the voltages of the 
transverse leads when /2 = -^3 = condition is imposed 
on their currents. 

The device Hamiltonian containing the Rashba SO 
term commutes with the unitary transformations which 
represent these symmetry operations in the Hilbert space 
'Ho ® 'Hs- (i) U{C2) exp {i^&z), which performs the 
transformation — > — (Tx, (Ty ~(^yy —* o'z and in- 
terchanges the leads 1 and 4 as well as the leads 2 and 3; 
(ii) U {(Jvx) ® exp (i^(Ty), which transforms the Pauli ma- 
trices Gx — > — CTx, Oy — > (Tj,, — —^z and interchanges 
leads 2 and 3; and (iii) U {(Tvy) <8) exp {i^&y) which trans- 
forms Gx ~* Gx, CTy ^CTj/, (^z ^ —(^z and exchanges 
lead 1 with lead 4. The Hamiltonian also commutes with 
the time-reversal operator if exp (i^Gy) [where K is the 
complex conjugation operator]. 

The effect of these symmetries on the spin-resolved 
charge conductance coefficients, and the corresponding 
spin conductances G^j^ , and Gfp expressed in terms 
of them through Eq. (jSJ, is as follows. The change in 
the sign of the spin operator means that spin-l becomes 
spin-| so that, e.g., G™ will be transformed into — G"p. 
Also, the time reversal implies changing the signs of all 
spin operators and all momenta so that G^^ = G~p 
is equivalent to G™g = —G°p^. Thus, invariance with re- 
spect to ?7((T„y)(g)exp {i^(Jy) yields the identities Gl^^'^ — 

G24 1 G^x^ ~ — ^24^' and G2\ — — G24. • These 
symmetries do not imply cancellation of G™^^ . How- 
ever, invariance with respect to U{ayx) ® exp {i^cfy) and 
IJ{C4) (g)exp (j jfTz) implies that Gjs^ = and Gj"'^ = 0. 

These symmetry imposed conditions simplify the gen- 
eral formula Eq. ^ for spin conductances of a perfectly 



clean Rashba SO coupled four-terminal bridge to 

Gf^ = 2G™'■^ + G™''^ (10a) 

GIh = G°^''^ (10b) 

GIh = G™*'^ (10c) 

where we employ the result V2/V1 = V3/V1 = 0.5 valid 
for a geometrically symmetric clean bridge. Because this 
solution for the transverse terminal voltages is violated 
in disordered bridges, its sample specific (for given im- 
purity configuration) spin conductances cannot be com- 
puted from simplified formulas Eq. (|1U|I . 

It insightful to apply the same symmetry analy- 
sis to the bridges with other types of SO couplings. 
For example, if the Rashba term in the Hamiltonian 
Eq. ^ is replaced by the linear Dresselhaus SO term 
f {Px^x — Py^y) due to bulk inversion asymmetryji no 
qualitative change in our analysis ensues since the two SO 
couplings can be transformed into each other by a unitary 
matrix {ctx + By) /\/2. In this case, the spin Hall response 
is signified by [/|]^ = — [/|]^ and [/|]'' = — [/|]^ compo- 
nents of the transverse spin current, while [/|]^ = [/|]^. 
For the Dresselhaus SO coupled bridge, the general ex- 
pression Eq. © simplifies to 

G^H = (11a) 
GIh - 2G™*•^-fG™''^ (lib) 
GIh = G°T'\ (11c) 

The qualitatively different situation emerges when both 
the Rashba and the linear Dresselhaus SO couplings be- 
come relevant in the central region of the bridge since in 
this case it is impossible to find spin rotation which, com- 
bined with the spatial symmetry, would keep the Hamil- 
tonian invariant while only transforming the signs of its 
spin matrices. Moreover, for such ballistic bridge the con- 
dition condition I2 =1^ = leads to V2/V1 = 1 — V3/V1 
solution for the voltages, whereas imposing the alterna- 
tive condition V2 = V3 generates non-zero charge currents 
fiowing through the transverse leads 2 and 3 together 
with the spin currents [for which no simple relations akin 
to Eq. © can be written in either of these cases] . 

III. TRANSVERSE PURE SPIN CURRENTS IN 
BALLISTIC BRIDGES 

Figures |21 El and 0] demonstrate that the spin Hall 
conductance is not universal — GsH{Ep,tsoT L,W) de- 
pends on the 2DEG parameters such as the density of 
charge carriers (i.e., the Fermi energy Ep), the Rashba 
SO coupling tso — ct/2a, and the system size L. Further- 
more, due to the sensitivity of spin dynamics in confined 
SO coupled ballistic systems to the boundaries and inter- 
faces,^"* it can also be affected by the measuring geometry. 

Nevertheless, we find in Fig. O that all square 2DEG 
samples of the size Lgo x LsOi where Lso is the spin pre- 
cession length, exhibit approximately the same G^j^ ~ 
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FIG. 2: (Color online) The Fermi energy dependence of the 
spin Hall G^Hy G^h a^nd spin polarization G|p conductances 
in ballistic four-probe bridges (Fig.0 where the central region 
is 2DEG of the size 60a x 60a. Within the 2DEG electrons 
are subjected to the Rashba SO coupling whose strength is 
tQo = a/2a = O.OUo. 



FIG. 3: (Color online) The dependence of the spin Hall G^^, 
G^H a-nd spin polarization G|p conductances on the Rashba 
SO coupling strength tso in the 2DEG of the size 60a x 60a 
attached to four ideal (with no spin tso = and charge inter- 
actions) leads. The unpolarized current injected through the 
longitudinal leads is composed of spin-^ and spin-| electrons 
at the Fermi energy Ef ~ — 3.8to. 



0.2e/47r. Therefore, we pay special attention to the in- 
tertwined effect of a and L brought about by the fact that 
SO couplings introduce a characteristic length scale into 
the bridge — the spin precession length Lso ~ 7r/2fcso 
over which spin precesses by an angle tt (i.e., the state 
It) evolves into |J,)). In the case of the Rashba SO cou- 
pling, 2fcso — 2m*a/h^ is the difference of Fermi wave 
vectors for the spin-split transverse energy subbands of a 
quantum wire. This quantity is the same for all subbands 
of the quantum wire in the case of parabolic energy- 
momentum dispersion so that single parameter 

^so = (12) 

characterizes the whole structure. For example, the 
mesoscale Lso sets the characteristic length for the evo- 
lution of the nonequilibrium spin polarization in the 
course of transport through the ballistiCf^^ as well as the 
diffusive SO coupled structures (which are sufficiently 
wide and weakly disordered^) where it plays the role 
of the disorder-independent D'yakonov-Perel' spin relax- 
ation lengthpiSiS^ 

In this section we study systematically the dependence 
of the spin Hall and spin polarization conductances on 
the first three basic parameters of the bridge, while leav- 
ing the effect of the disorder strength W for Sec. IIVI (cf. 



Fig. O , and the influence of the measuring geometry de- 
termined by the attached leads to Sec.El(cf. Fig. EJ. 

A. Fermi energy dependence of the spin Hall 
conductance 

The spin conductances plotted in Fig. El are an odd 
function of the Fermi energy and, therefore, have to van- 
ish GIh{Ef = 0) = 0, GVp{EF = 0) = at the half-filled 
band Ep — 0. This feature is a consequence of the spin 
current being defined as the difference of charge currents 
of spin- 1 and spin- J, electrons and particle-hole symmetry 
of the tight-binding Hamiltonian Eq. ||SJ). That is, the 
spin current carried by electrons above the half-filling 
Ep > can be interpreted as the propagation of posi- 
tively charged holes which move in the opposite direction, 
and have opposite spin, to that of electrons. To evade 
artifacts of the tight-binding energy-momentum disper- 
sion, which can enhance the spin Hall conductance as 
we approach the band center, we highlight in Fig. H the 
values of spin conductances at the Fermi energy Ep — 
—3.8to chosen for our subsequent analysis. When zero- 
temperature unpolarized charge quantum transport is de- 
termined by the states at this Ep, which is close to the 
bottom of the band Eh = — 4to, the injected quasiparti- 
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FIG. 4: (Color online) The finite-size scaling of the spin Hall G^^, G^^ and spin polarization Ggp conductances in ballistic 
four-probe bridges with the Rashba SO coupling strength: (a) tso ~ O.Olto; and (b) tso = O.lto. The corresponding spin 
precession length denoted on the graphs is: (a) Lso ~ 157a; and (b) Lso ~ 15.7a. 



cles have quadratic and isotropic energy-momentum dis- 
persion which characterizes the Hamiltonians in effective 
mass approximation, such as the Rashba one in Eq. 
Also the Fermi wavelength corresponding to the Fermi 
energy measured from the band bottom Ep — Ef, = 0.2to 
is much greater than the lattice spacing a so that possible 
artifacts of the discretization are avoided. 



B. Rashba SO coupling dependence of the spin 
Hall conductances 



Considerable interest for exploiting the Rashba SO 
coupling^j^ for semiconductor spintronics apphcations 
stems from the possibility to tune its strength via an ex- 
ternal gate electrodej^ thereby manipulating spin solely 
by electrical means. A surprising early result in the 
theory of the intrinsic spin Hall effect is apparent 'uni- 
versality' of = e/SiT (obtained from the linear re- 
sponse theory for the clean Rashba Hamiltonian of an 
infinite 2DEG) in the sense that it does not depend on 
the strength of the SO couplingi^ However, when scat- 
tering of impurities is taken into account in the limit 
a — !■ 0, one recoverg^^> the physically expected result 
lima^o cTsH 0. 

Our exact treatment of transport of non-interacting 
quasiparticles through a clean finite-size system does not 
face any technical impediments in locating the lower limit 
on the strength of SO coupling capable of inducing the 



non-negligible spin Hall conductance G^j^, as shown in 
Fig. 121 Although realistic Rashba SO coupling strengths 
in current experiments can be tuned within the range 
O.Olio ^ ^so S 0-Uo (see Seem, '^^ also show in the 
same Figure the upper limit of large values of tso beyond 
which all three components of the transverse spin current 
/| vanish due to the carrier reflect ionSLS^ at the interface 
between the ideal lead with ^so = and the sample with 
strong SO coupling tso ~ ^o- The maximum G^^(tso) 
for square shaped 2DEGs attached to ideal semi-infinite 
leads is obtained for samples of the size Lso x -^SO, which 
connects this dependence GJjy(iso) with the finite-size 
scaling properties of Glu{L) discussed in the next sec- 
tion. 



C. Finite-size scaling of the spin Hall conductances 

FigureElemphasizes the importance of the spin preces- 
sion length scale for the spin Hall effect in ballistic semi- 
conductor nanostructures. That is, the spin Hall con- 
ductance GJ^(L) is increasing non-monotonically with 
the system size for L < LsO; attaining the maximum 
value GfjLf(Lso) — 0.2e/47r. In such L x L samples with 
L < Lso, the other spin Hall conductance Ggfj{L) is 
negligible. The emergence of [/|]^ component of the spin 
Hall current can be understood heuristically by invok- 
ing the semiclassical picture involving the transverse SO 
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force operator— 

Fso ^ ^3 (p X z) d^'^' ^ ^ 

and studying its effect on the propagation of spin- 
polarized wave packetSii^ This spin-dependent terms 
here are generated by the Rashba SO couphng term in the 
single-particle Hamiltonian of a clean 2DEG. Its expecta- 
tion values in the wave packets j'l') (X) jcr) (spin-polarized 
along the z-axis, (7z\^) = +|T) and az\l) = shows 
that spin-t and spin- J, injected electrons will be deflected 
in opposite transverse directions (e.g., spin-| is initiallyi^ 
deflected to the right). Moreover, their spins are forced 
into precession since injected |t), |i) states are not the 
eigenstates of the Zeeman term & ■ Bfl,(p), where the 
Rashba effective magnetic field 'Bj^{p) remains nearly 
parallel to the y-axis due to the transverse confining po- 
teiitia.&^Vconiiy). 

The precession of the defiected spins is responsible for 
the oscillatory character— of the SO "force" (viewed as 
the expectation value of the SO force operator in the 
spin-polarized wave packet states) which changes sign 
along the wire. Thus, such a^-dependent spin-deflecting 
"force" can induce the change in the sign of the spin 
Hall current as a function of the system size, as shown 
in Fig. Efb) for strong SO coupling tso — O.lto. In the 
case of weaker SO couplings tso ^5 0.04io, the spin Hall 
conductance oscillates for L < Lso while remaining 
positive, as exemplified by Fig. |3|[a). Note that in the 
convention of multiprobe spin current formulas Eq. Q), 
positive spin Hall conductance means that spin current 
/| is flowing out of the lead 2 because spin-t electrons 
are deflected to the right (i.e., toward the electrode 3 if 
injected from the electrode 1) and spin-| electrons are 
deflected to the left, as expected from Eq. 

The semiclassical picture of the SO "force" also ex- 
plains the symmetry properties of the spin Hall current 
[/|]^(— y) = [/|]^(F) with respect to the voltage bias 
reversal [i.e., the reversal of the momentum in Eq. 
or the reversal of the sign of the Rashba SO coupling, 
[^2]^(~o;) = [^l]^!*^)- These two features make it possible 
to differentiateii*i^ between the z- and the x-component 
of the spin Hall current in the transverse leads since 

mn-v) = mnv) and mn-a) = -mna), which 

stem from the properties of the effective magnetic field 
B_r(p) inducing the spin precession under these transfor- 
mations. 

The y-component of the transverse spin current satis- 
fies [/|]'' = [/|]'', signaling that Gf^ is of completely dif- 
ferent origin. It stems from {Sy(r)) ^ nonequilibrium 
spin accumulation, which has the same sign on the lat- 
eral edges of 2DEG,^^ induced when unpolarizcd charge 
current is flowing through the Rashba SO coupled 2DEG 
attached to two electrodesiii*22iS& Thus, when transverse 
leads are connected to the lateral edges the 2DEG, the 
spin-dependent chemical potentials^ on the edges will 
push the spin current [/|]*' into the leads. Since this 
chemical potential is the same on both lateral edges, con- 



necting the edges by a wire would not lead to any net 
spin flux through it, in sharp contrast to the currents 
[I2Y and [J|]^ that will transport spin through such wire 
as a signature of the spin Hall effect (s) phenomenology. 



IV. TRANSVERSE PURE SPIN CURRENTS IN 
DISORDERED BRIDGES 

The most conspicuous difference between the "old" ex- 
trinsic^iSi effect (in paramagnetic metals or semiconduc- 
tors without SO splitting of quasiparticle energies) and 
the "new" intrinsic spin Hall effect (in semiconductors 
with sufficiently large SO splitting of quasiparticle en- 
ergies) is that the former vanishes in the ballistic limit, 
while the later persists even when no skew-scattering at 
impurities takes place. ^'^ However, the distinction be- 
tween the two spin Hall effects turns out to be ambigu- 
Qyg27j28 .y^rj^gjj so energy splitting akp [kp is the Fermi 
wave vector) is smaller than the disorder induced broad- 
ening of the energy levels h/r (where r is the transport 
lifetime). 

Moreover, the resilience of the intrinsic effect to scat- 
tering off static impurities has become a major issue in 
current debates over the observability of spin Hall cur- 
rent in realistic samples of 2DEG with the Rashba SO 
interaction2Si2LS2422i2a or in hole-doped bulk 3D semi- 
conductors >2L2i Early (lowest order) perturbative treat- 
ment of the semiclassical Boltzmann diffusive transport 
in infinite homogeneous Rashba spin-split 2DEG has re- 
vealed that UsH could survive disorder effects with the 
proviso that akpr/h 3> 1, while being gradually reduced 
from the 'universal' value e/87r with increasing disor- 
der strength— However, when vertex corrections are in- 
cluded in the perturbative expansion, the intrinsic effect 
turns out to be suppressed asH — > at arbitrarily small 
disorder akpr/h ^ 00, as revealed by a multitude of 
transport approaches'^ applied to any system with linear 
in momentum SO energy splitting (note that such can- 
cellation induced by the ladder vertex corrections does 
not occur when SO coupling contains higher-order mo- 
mentum termSi^). 

The vanishing intrinsic spin Hall current density jy — > 
is also found at the weak localization level in the pertur- 
bative expansion in small parameter— l/kp£, as well as 
in the bulk of 2DEG (infinite in the transverse direction) 
attached to two massive electrodes in the longitudinal di- 
rection where, nevertheless, macroscopic inhomogeneities 
can induce non-zero jy ^ within the spin relaxation 
length Lso wide region around the electrode-2DEG in- 
terfaces— Finally, recent reexamination^^ of these re- 
sults in a vast range of ratios of the SO coupling induced 
akp and the disorder induced h/r energy scales confirms 
that asH — > is indeed suppressed in both the clean 
akpr/h — > 00 and dirty limits akpr/h — > 0, while being 
able to attain a non-zero optimum value at akpr/h ~ 10. 

Here we shed light on the effect of disorder on the 
spin Hall current /|, whose maximum value is set by 
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FIG. 5: (Color online) The effect of the spin-independent 
static disorder on the spin Hall {G^u), (G^jf) and spin po- 
larization (Gsp) conductances in the four-probe mesoscopic 
bridges with different strengths of the Rashba coupling tso- 
In the weak disorder limit, semiclassical mean free path is 
£ ^ 21.5atl/W^ , while in the shaded range of W both the 
disorder-averaged (G) and the typical e^'"*^^ two-probe charge 
conductance of the SO coupled 2DEG (attached to two leads) 
diminish below ~ 0.1 (2e^//i) due to strong localization effects. 



the SO coupling effects in ballistic transport regime 
through multiprobe mesoscopic structures, by intro- 
ducing spin-independent static random potential Em S 
[-W/2, W/2] into the 2DEG region of Fig. [Hand study- 
ing the disorder-averaged conductances (Gf^), {G^,p), 
(GJ^) in the crossover from the quasiballistic to the lo- 
calized transport regime. The inclusion of all transport 
regimes is made possible by employing the exact single- 
particle spin-dependent Green function Eq. |(7J) which en- 
compasses all quantum-interference effects at arbitrary 
W and ^so in a finite-size device, rather than treat- 
ing only the lowest order (semiclassical) effects of the 
disorderiSSiSLSSiSS or weak localization quantum correc- 
tions.^*^ Figure |S1 suggests that mesoscopic spin Hall con- 
ductances are unaffected by weak disorder, and they 
will gradually diminish toward negligible values only at 
akph/r = (tso^) / (toO.) — 0.1 which is deep inside the 
diffusive metallic £ < L <^ regime. Here the semiclassi- 
cal mean free path £ w 21.5at^/W^ (at the selected Fermi 
energy Ep — — 3.8to) is applicable within the Boltzmann 
transport regime^^ £ > a. 

The often quoted mantra — all quantum states of dis- 
ordered non-interacting electrons in two-dimensions are 
localized — means in practice that the localization length 



f < 00 is finite at arbitrary disorder strengthn^ Thus, 
the two-probe (charge) conductance of a sufficiently large 
L ^ ^ 2DEG with impurities will decay exponentially 
fast G ~ The two exceptions are quantum Hall 

2DEG (where delocalized states exist in the center of a 
Landau leveliS) and 2D systems where strong enough SO 
coupling can induce the metallic phase ^ — > oo at weak 
disorder. '^'^ We delineate in Fig. [31 the boundaries2S of 
the localization-delocalization transition to demonstrate 
that spin Hall conductances will vanish upon increasing 
disorder before 2DEG is pushed into the realm of strong 
localization where GsH — > is trivially expected. 



V. MESOSCOPIC SPIN HALL CONDUCTANCE 
VS. INTRINSIC SPIN HALL CONDUCTIVITY 

Both the intrinsic and the mesoscopic spin Hall ef- 
fect originate from the SO coupling terms in Hamilto- 
nians of clean semiconductor systems. Since conductiv- 
ity and conductance have the same unit in two dimen- 
sions, one might naively expect that G^^ = OsH since it 
is tempting to connect the total spin Hall current with 
the spin current density integrated over the cross sec- 



tion, /| 



3lL. 



and use = {Vi — V^j/L to find that 



G^jLf = UyL) I {ExL) = ffsH should not scale with the sys- 
tem size. If, on the other hand, mesoscopic spin Hall cur- 
rent in Rashba SO coupled devices is due to the edge spin 
currents near the contacts with the longitudinal leadS)S& 
the same arguments would lead to G^fj (x Lgo/L which 
decreases with the 2DEG size. 

In contrast to these naive expectations, Figure |21 re- 
veals more complicated scaling behavior of Gljj{L) in 
the "mesoscopic" regime where spin Hall conductance 
oscillates for L < isO: and reaches asymptotic value 
Glu ~ 0.1e/47r (up to small oscillations around it due 
to phase-coherent and ballistic nature of transport) in 
the "macroscopic" regime L 3> ^so- The semiclassical 
picture of the deflection of spin densities involving the SO 
"force" Eq. (|13l) offers an explanation of these oscillations 
of the spin Hall conductance as being due to the change 
in sign of such "force" deflecting the spin which at the 
same time is precessing, as elaborated in Sec. IIII CI The 
decay of the SO "force" magnitudei^ while the spin is 
moving along the wire, which arises due to spin decoher- 
ence in ballistic SO coupled systemS)22iS4 is responsible 
for the saturation of spin Hall conductance in large sam- 
ples with strong SO coupling. 

A closer look reveals further fundamental differences 
between G^^ and asH- The intrinsic spin Hall effect 
which is driven by an external electric field penetrating 
infinite homogeneous systems in the clean limit is 
essentially a semiclassical phenomenon where spin cur- 
rent is generated by the anomalous velocity (due to the 
Berry phase in momentum space) of Bloch wave packets 
without requiring the shift of the electron distribution 
function from equilibrium.'^'^ Such unusual properties of 
jy carried by the whole SO coupled Fermi sea, which 
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depends only on the equilibrium distribution function 
and spin-split band structure)i^ have lead to arguments 
that the intrinsic spin Hall current is an equilibrium cur- 
rent which does not actually transport spin between two 
points in space.'^^'^'^ In fact, such jy ^ 0, which does not 
induce spin accumulation or can be employed for spin in- 
jection, is found in SO coupled systems without any ap- 
plied electric field'^^ or obvious sources of dissipation^ 
which is compatible with time-reversal invariance since 
jy does not change sign under the time-inversion t —t 
transformation. 

On the other hand, the mesoscopic spin Hall current 
/| is a Fermi-surface quantity at zero temperature T — > 
[i.e., the contribution to GsH in Eq. |(SJ) from Green func- 
tions Eq. lO evaluated at energies E < Ep is zero] and 
a genuine nonequilibrium response because no total spin 
currents can flow throughout the leads of a multitermi- 



nal device in equilibriui 



,19.20 



const, (as discussed 
in Sec. Moreover, the induction of J| can never be 
"dissipationless"i^ since GsH in Eq. ||SJ) is expressed 
in terms of the spin-resolved charge conductance coef- 
ficients whose non-zero values, even in perfectly clean 
systems where G^^ > is set entirely by the sample 
geometrjii^ and interfaces at which SO coupling changes 
abruptly^SiiSS encode the information about dissipation 
occuring in remote huge reservoirs thermalizing electrons 
to ensure the steady-state transport. That is, in the Lan- 
dauer setup current is limited by quantum transmission 
through a potential profile while power is dissipated non- 
locally in the reservoirsji^ Note also that in the four- 
terminal devices of Fig. ^ the external bias voltage only 
shifts the relative chemical potentials of the reservoirs 
into which the longitudinal leads eventually terminate so 
that electrons do not feel any electric field in the course of 
ballistic propagation through clean 2DEG central region. 

The inability to connect bulk conductivity (which re- 
lates local current density to the electric field, j — crE) to 
conductance measured in experiments (which relates the 
total current to the voltage drop, / = GV) is encountered 
in some charge transport situations as well: (i) in the 
ballistic regime £ ^ L conductivity a does not exist as a 
local quantity and only the conductance G plays a role; 
(ii) in phase-coherent diffusive conductors quantum cor- 
rections to the conductivity emerge which are non-local 
on the dephasing scale that can be much greater than 
the mean free path (typically < 1/im below which the 
inelastic processes become suppressed at low enough tem- 
peratures T <^ IK), so that concept of local conductivity 
loses its meaning. Although our devices in Sec. IHII are 
both in the ballistic £ ^ L and phase-coherent L < 
regime, the principal obstacle in connecting the intrinsic 
spin Hall conductivity UsH and mesoscopic spin Hall con- 
ductance GIjj lies in the fact that spin is not conserved 
in SO coupled systems. Thus, the plausibly defined 
spin current-density operator = j{o'iVk + Vk^i) (sym- 
metrized product of the Pauli matrices and the velocity 
operator which yields a Hermitian operator), employed 
in different computational schemes for the bulk spin 
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FIG. 6: (Color online) The spin Hall conductance Ggu of 
a 2DEG of size 15a x 15a with the Rashba SO coupling 
(dark zone), which is attached to four ideal leads with no 
Rashba coupling (solid line) or four leads where SO coupling 
is switched on adiabatically via a linear function within a 
finite region (light zone) of length 15a of the leads (dotted 
line). Panel (a) compares the Fermi energy dependence of 
Glif for these two difi'erent measuring setups, while panel (b) 
compares their SO coupling dependence. 



Hall conductivity; '^'i^i^^i^'^'i^^i^^i^"'^^ lacks rigorous theoret- 
ical justification and standard physical interpretation be- 
cause it does not satisfy the continuity equation with the 
spin density operator jiSL^S, 

While the relation of the spin current density jy (as 
the expectation value of the corresponding "controver- 
sial" spin current-density operator) flowing through the 
SO coupled system to real spin transport and spin ac- 
cumulation is far from obvious, the pure spin current 
^(/J — I2), which we define within the asymptotic re- 
gion of the leads with a = 0, is conserved quantity so 
that /| = const, does not change on different transverse 
cross sections throughout the lead 2. Such pure spin cur- 
rents flowing through the region with no SO coupling 
have transparent physical interpretation: If all spin-^ 
electrons move in one direction, while an equal number 
of spin-| move in the opposite direction, the net charge 
current vanishes while spin current can be non-zero. In 
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fact, they have been created and detected in recent opti- 
cal pump-probe experimentsjS 

Even in the semiclassical transport regimes at higher 
temperatures (where L > L^), the presence of SO cou- 
phng emphasizes the demand to treat the whole device 
geometry when studying the dynamics of transported 
spin densities. For example, the decay of nonequihb- 
rium spin polarizaions in baUistic or disordered quantum 
wires is highly dependent on the transverse confinement 
effectsS^ or chaotic vs. regular boundaries of clean quan- 
tum dotsiSi This is to be contrasted with the conven- 
tional D'yakonov-Perel' spin relaxation mechanismJ- in 
unbounded diffusive systems where the decay of spin po- 
larization is determined solely by the SO coupling and 
elastic spin-independent scattering of charges on the im- 
purities in the bulkJi Also, the eigenstates^^ of SO cou- 
pled wires substantially differ from the ones of the infi- 
nite 2DEG since Bfl;(p) is almost parallel to the trans- 
verse directioniSiSiiSS (in contrast to the infinite 2DEG 
where no unique spin quantization axis exists^). These 
transverse confinement effects are found in Sec. 1111 CI to 
be responsible for the non-zero [/|]^ ^ component of 
the mesoscopic spin Hall current, which is quite different 
from jy = property of the intrinsic spin Hall effect in 
infinite homogeneous systems. 

The mesoscopic transport techniques, developed to 
treat the whole measuring geometry as demanded by 
quantum coherence effects and non-local nature of trans- 
port measurements in phase-coherent devices, are well- 
suited to handle all relevant details of the spin Hall 
bridges, as applied in Sec. ^ to spin transport. We in- 
vestigate this issues further in Fig. by studying the ef- 
fect of measuring geometry on the maximum value of the 
spin Hall conductance in the four-probe bridges with 
Lso X -^so 2DEG, where we employ the leads contain- 
ing finite region within which the Rashba SO coupling 
is switched on adiabatically (via linear function) from 
^SO = to the value it acquires in the 2DEG. Compared 
with our standard setup from Fig. the usage of such 
electrodes would enhance the spin Hall effect since reflec- 
tion at interfaces where iso changes abruptly is greatly 
reduced. 



VI. CONCLUDING REMARKS 

In conclusion, we have delineated features of a novel 
type of spin Hall effect in four-terminal mesoscopic struc- 



tures where unpolarized charge current driven through 
the longitudinal ideal (with no spin and charge interac- 
tion) leads attached to clean semiconductor region with 
strong enough homogeneous SO coupling induces a pure 
spin current in the transverse voltage probes (with zero 
net charge flow through them). The spin carried by the 
transverse spin Hall current in devices where electrons 
are subjected to the Rashba-type of SO coupling within 
the central region has both out-of-plane and in-plane non- 
zero components of its polarization vector. The spin Hall 
current depends on the strength of the Rashba SO cou- 
pling. Furthermore, the maximum value of the spin Hall 
conductance Gf^ is always achieved for the square sam- 
ples of the size Lso x Lgo, where spin precession length 
Lso can be tuned by changing the Rashba SO coupling 
via the gate electrode covering the 2DEG. 

Although apparently similar to recently predicted in- 
trinsic spin Hall effect (as a semiclassical phenomenon 
in infinite homogeneous clean SO coupled systems), the 
mesoscopic pure spin Hall current predicted here has fun- 
damentally different properties: It is a genuine nonequi- 
librium and Fermi-surface quantity which depends on the 
whole device (i.e., interfaces and boundaries) and mea- 
suring geometry (e.g., it can be enhanced by the leads 
where SO coupling is adiabatically switched on within a 
finite region). Its non-trivial finite-size scaling regimes 
in samples smaller and larger than Lgo highlight the im- 
portance of processesi^ involving spin dynamics on this 
mesoscale for the generation of spin Hall current in bal- 
listic multiprobe nanostructures. In contrast to the bulk 
intrinsic spin Hall conductivity of the Rashba spin-split 
2DEG, mesoscopic spin Hall conductances are able to 
survive weak scattering off impurities, gradually decay- 
ing within the metallic diffusive transport regime. 

Note added. — Upon completion of this work we have 
become aware of Ref. [s^ where similar mesoscopic ap- 
proach to spin Hall effect has been undertaken for dif- 
ferent four-probe structure — an infinite Rashba SO cou- 
pled wire with two transverse ideal leads attached — and 
analogous conclusions have been reached regarding its 
resilience to disorder effects. 
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